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Abstract 

The fermion mass matrix, in addition to having eigenvalues (masses) 
which run, also changes its orientation (rotates) with changing energy 
scales. This means that its eigenstates at one scale will no longer be 
eigenstates at another scale, leading to effects where fermions of dif- 
ferent flavours can "transmute" into one another. In this paper, the 
implications of a rotating mass matrix are analysed and possible trans- 
muation effects are investigated both in the Standard Model (SM) and 
in the so-called Dualized Standard Model (DSM) that we advocate, 
arriving at the conclusion that some transmutational decays such as 
ip — > ut, T — ► ut or 7T — ► e/i may be within experimental range, 
if not immediately, then in the near future. 



1 Introduction - the Rotating Mass Matrix 



By a rotating mass matrix, we mean one which undergoes unitary transfor- 
mations through scale changes. 

As most quantities in quantum field theory vary under changing scales, 
in particular the masses of particles or in other words the eigenvalues of mass 
matrices, it is no surprise that the orientations of mass matrices may also 
change. 

Indeed, even in the Standard Model as usually conceived, nontrivial mix- 
ing between up and down fermion states will generally induce rotating mass 
matrices. That this is so can be seen from the following renormalization 
group equations: 

16vr 2 ^ = l(Urf - DD^U + (S u - A U )U, (1.1) 

16vr 2 ^ = |pl>t _ uu\ )D + ( Ed _ A d )D. (1.2) 

satisfied by the mass matrices U and D of respectively the up and down 
fermions Nontrivial mixing between up- and down-states means that 
the matrices U and D are related by a nondiagonal matrix D = VUV^ 
or U — V^DV, so that U and D cannot be simultaneously digonalized. 
Suppose now at some t we diagonalize U in ( |1.1|) by an appropriate unitary 
transformation. All other terms in (|1.1|) are then diagonal but not the term 
DD^U. Hence, this term will necessarily de- diagonalize the matrix U on 
running by ( |1 . 1\) to a different t, or in other words the matrix U will rotate 
with changing t 0. Similar arguments hold also for the rotation of the 
matrix D. Given that quarks have long been known to have a nontrivial 
(CKM) mixing matrix ||, and that recent experiments on neutrino oscillation 
||, [|, Q strongly suggest a nontrivial (MNS) mixing ]7) also for leptons, we 
have to conclude that the mass matrices for both quarks and leptons will 
rotate with changing scales even if there are no forces and interactions in 
nature other than those currently studied in the Standard Model. 

However, there is clearly a possibility — perhaps one might even say a 
strong theoretical reason to suppose — that forces may exist in nature other 
than those now conventionally studied in the Standard Model, which can give 
further rotations to the fermion mass matrices. Indeed, in the usual formula- 
tion of the Standard Model, the fact that there are 3 generations of fermions 
and that they mix is taken as an input from experiment, with the result that 
the framework depends on a large number of empirical parameters, about 
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three-quarters of which are traceable to the mystery of generations. One 
could thus hope that in future when the generation puzzle is solved, some 
or even most of the empirical parameters appearing in the present Standard 
Model will be predictable. One favourite and perhaps most natural assump- 
tion for theoretical attempts in this direction is that generations originate 
as a broken 'horizontal' gauge symmetry ||. If that is the case, then the 
interactions associated with this symmetry which mix the generation index 
are likely also to rotate the fermion mass matrices in generation space. This 
rotation would be over and above that driven by the nondiagonal CKM or 
MNS matrices via the mechanism described in the preceding paragraph. 

Whatever its origin, however, the rotation of fermion mass matrices is 
theoretically of the same standing as the running of masses and coupling 
constants, and if testable by experiment would lend equal support to the 
quantum field theory framework. Moreover, if experiment shows that the ro- 
tation is different from simply that given by renormalization group equations 
fll.lj) and ( |1.2|) via the empirical mixing matrix, then it suggests that there 



are forces at work other than those currently studied in the Standard Model 
which are of the type considered in the preceding paragraph. This means 
that we would have opened up a new window for investigations furthering 
our long-standing quest for the origin of fermion generations. And since, as 
we shall see, mass matrix rotation leads to an entirely new and very distinc- 
tive category of physical phenomena, there is a fair chance of its effects being 
observed by experiment in the not too distant future. The purpose of the 
present paper is to explore such possibilities. 



2 Specification of Fermion States 

The fact that the mass matrix rotates poses immediately a question of physi- 
cal interpetation at a rather basic level, namely the question of how to define 
the state vectors of the various fermion states, and hence the mixing ma- 
trices between them. For the familiar case of a non-rotating mass matrix, 
one defines the state vectors of the 3 generations as its eigenvectors, which 
are by assumption scale-independent and, since the matrix is hermitian, are 
also mutually orthogonal as they should be if they are to represent indepen- 
dent physical entities. The mixing matrix between the up and down fermion 
types, being the overlap matrix between 2 scale-independent orthonormal 
triads of state vectors, one for each type, is then also scale-independent and 
automatically unitary. Only the eigenvalues of a mass matrix then run, and 
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the masses of the 3 states (generations) can be defined as the running values 
each evaluated at the scale equal to its value. However, when the mass ma- 
trix rotates, although one can still diagonalize the matrix at any scale by a 
triad of mutually orthogonal eigenvectors, this triad will be scale-dependent, 
and it is not immediately clear which vectors at which scale(s) are to be 
identified as the state vectors of the 3 generations. Although the ambiguity 
may not in all cases be numerically significant in view of present limitations 
in experimental accuracy, it still has to be resolved as a matter of principle 
against the day when better accuracy is achieved. 

Now, the state vector of a physical state like the muon is normally, we 
think, taken as a scale-independent concept, or otherwise we are likely to 
meet with some awkwardness in its physical interpretation. For example, 
suppose we were to define the muon state as the eigenvector with the second 
highest eigenvalue of the charged lepton mass matrix at any scale. Then since 
the matrix rotates, the [i state vector will point at different scales in different 
directions in generation space. In that case, the muons obtained say from an 
energetic beam of charged 7r's, which we know decay almost entirely into /i's 
at rest, will no longer appear as purely muons when it hits the target, but as 
a linear combination of e, fi and r, which is at variance with what is usually 
understood. Although one can in principle insist on defining fermion states 
as the eigenstates of the rotating mass matrix at every scale, and so long as 
one is consistent in their interpretation one would arrive in the end at the 
same physical results, such a procedure would seem to be rather inconvenient. 
In this paper therefore, we opt for a scale-independent definition of fermion 
state vectors, in which case it would be incumbent upon us to specify exactly 
the scale at which each of these vectors are to be defined. 

If the rotation of the mass matrix is negligible to a certain approximation 
below some scale, then one may define to this approximation the state vectors 
as the orthonormal triad of eigenvetors of the mass matrix taken at that scale. 
This, though seldom stated explicitly is, we think, the tacit criterion usually 
adopted in the literature which we shall refer to in this paper as Fixed Scale 
Diagonalization FSD. It is, of course, at best an approximation for the mass 
matrix cannot stop rotating abruptly, and that approximation, as we shall 
see later, may not be all that good in certain circumstances. Besides, this 
criterion involves a degree of arbitrariness in choosing a certain fixed scale 
to effect the diagonalisation of the mass matrix, and begs the question of 
principle why it should be that particular scale and not some other. 

In view of this, one may be tempted instead to do the following. As one 
defines the mass of a state as the eigenvalue evaluated at the scale equal 
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to its value, one may try to define as its state vector the corresponding 
eigenvector also at the same scale. This would be more democratic and less 
arbitray than the FSD criterion but it will not work. Consider for example 
the [/-type quarks. Following the above proposal, one would then define the 
t state vector as the eigenvector with the largest eigenvalue 777,3 of the matrix 
evaluated at the scale /x = m 3 (/i), the c state vector as the eigenvector 
with the second largest eigenvalue m% at the scale /i = m 2 (/i), and the u 
state vector as the eigenvector with the smallest eigenvalue mi at the scale 
/j = mi(/i). These 3 state vectors, however, will not be mutually orthogonal, 
for although the 3 eigenvectors of the hermitian mass matrix are mutually 
orthogonal when all evaluated at the same scale, they have no reason to 
be so when evaluated each at a vastly different scale since the matrix by 
assumption rotates with changing scales. It thus contradicts the assertion 
that they represent 3 independent physical entities; in particular, it would 
imply that the mixing (CKM or MNS) matrix, being a transformation matrix 
between the 2 triads of up and down state vectors, is not unitary as it ought 
to be. 

There is, however, a working criterion for defining state vectors each at 
its own mass scale along the lines suggested in the above paragraph which 
takes account of the rotation between the different mass scales yet still gives 
mutually orthogonal state vectors and hence a unitary mixing matrix. This 
was first proposed by us |9j in connection with a scheme for quark and lepton 
mixing which we call the Dualized Standard Model (DSM), but the criterion 
can in fact be applied to define state vectors in any scheme with a rotating 
mass matrix HlO| . For 3 generations, the criterion goes as follow. We run 
the mass matrix m down in scale until we have for its highest eigenvalue 7773 
a solution to the equation \i = m 3 (/i). This value at this scale we define 
as the mass 777.3, an d the corresponding eigenvector the state vector v 3 of 
the heaviest generation. Below that energy, the state 3 no longer exists as 
a physical state, and only the two lighter generations survive, whose state 
vectors have to be orthogonal to V3. We define then the mass matrix at 
energies below 7773 as the 2x2 submatrix rh of m in the subspace orthogonal 
to V3. To find now the mass and state vector for generation 2, we follow with 
m the same procedure used above with m for generation 3 and run m down 
in scale until we find a solution to the equation /i = 7772 (/i), which value we 
call the mass 777 2 and the corresponding eigenvector at that scale the state 
vector v 2 of the generation 2. The state vector of the lightest generation 1 is 
now also defined as the vector orthogonal to both v 3 and v 2 , while the mass 
of 1 is obtained by repeating the same procedure, namely by running down in 
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scale the remaining lxl submatrix, namely the expectation value (vi|ra|vi), 
until its value equals the scale. In this way, each mass is evaluated at its own 
appropriate scale while the physical state vectors of the 3 generations are 
mutually orthogonal as they should be. Furthermore, the mixing matrix, 
taken as the overlap matrix between the triads of state vectors so defined 
for the up and down states, will be both unitary and scale- independent. 
This criterion for defining state vectors, which is the only one we are aware 
of with the required properties, we shall refer to in future as Step-by-Step 
Diagonalization SSD. 

3 Fermion Transmutation as Consequence 

Whether one accepts the SSD proposal for defining physical fermion state 
vectors or chooses to ignore these subtleties and adopts the approximate 
criterion of FSD, the fact remains that the mass matrix will continue to 
rotate at high energy, and so long as the state vectors are defined as the 
eigenvectors diagonalizing the mass matrix at some prescribed scale(s), they 
will generally no longer diagonalize the mass matrix at some higher energy 
scales. Now, the fermion mass matrix appears for example in the fermion 
propagator of Feynman diagrams, and if it is nondiagonal, so also will be the 
transition amplitudes represented by these diagrams. Physically, this means 
that transitions will be induced in which a particle of one generation converts 
into a particle of another generation. For example, the diagram of Figure [I], 




+ 
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Figure 1: Transmutation of e into \i in e + e collision. 
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where a tilde denotes the occurrence of a nondiagonal fermion mass matrix, 
will have an off-diagonal element linking the e state vector to the /i state 
vector, leading to a e + e~ — > n + e~ transition. 

Such a conversion of e to fi is quite distinctive and singular, differing from 
that induced, for example, by a flavour- changing neutral current (FCNC) 
exchanged between the initial e + and e~ . The latter will normally produce 
in addition to the fi + at the upper vertex also a fi~ at the lower vertex, 
giving rise to a "double conversion", although under certain circumstances 
where the exchanged FCNC bosons themselves mix, "single conversion" can 
also occur, as for example in the situation considered in [13|. However, even 
in that last case, the conversion rate depends on the masses and couplings 
of the exchanged FCNC bosons which are still rather elusive quantitites, 
whereas the conversion effect considered here is due only to the rotation of 
the mass matrix, and may therefore be quite restricted. We propose therefore 
to call the latter "transmutation" to distinguish it from other conversion 
phenomena. 

The effect of transmutation is so unique as a consequence of the rotating 
mass matrix that if it occurs with appreciable cross sections, it should be 
identifiable in experiment with little difficulty At first sight, therefore, the 
prediction looks alarming, for the conservation of such quantities as muon 
number have been checked already to high accuracy, and their wholesale 
violation as suggested above seems hardly likely to survive the existing ex- 
perimental bounds. On closer examination, however, the conclusion becomes 
less obvious, mainly for the following reason. At low energy, near the scale (s) 
at which the mass matrix is diagonalized to define the physical fermion states, 
as explained in the preceding section, the deviations of the mass matrix from 
diagonality are small, and so also will be the cross sections for transmuta- 
tion. At high energy, on the other hand, where the off-diagonal mass matrix 
elements can become sizeable, it turns out that transmutation cross sections 
will be suppressed, as can be seen in the following example. 

One of the simplest transmutation processes that one can think of is 
what one can call the photo-transmutation of leptons, which to leading order 
is given just by the Feynman diagrams in Figure |2] which are the same as 
for ordinary Compton scattering except that the fermion line here carries 
a rotating and generally non-diagonal mass matrix. Consider first the s- 
channel diagram (a) in which the lepton propagator + i/k — m) appears, 
carrying the non-diagonal mass matrix m where transmutation originates. 
For s = (p + k) 2 large, we can expand the propagator in powers of m/y/s, 
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Figure 2: Feynman diagrams for photo-transmutation of leptons. 
obtaining symbolically: 

1 1 , m . . 

l + -7r, 3.1 



where one sees that the off-diagonal terms leading to transmutation is sup- 
pressed by a factor of 1/ ^fs. Similarly, the lepton propagator in the diagram 
(b) for u large can be approximated as: 

1 1 , m „ , 

l + -p=f), (3-2) 



foi - - m) (/* - f) 

leading to transmutation effects suppressed by a factor 1/ J\u\. In the region 
where \u\ is small, off-diagonal contributions can remain sizeable, but the 
size of this region itself decreases as 1/s so that the integrated effect remains 
small. 

For these reasons, transmutation effects arising from the rotating mass 
matrix are seldom very large and the question whether they can already be 
ruled out by existing data, and if not, how soon, and under what conditions, 
they will be verifiable, is not immediately obvious and can only be answered 
by detailed analysis. The answers will depend on the process investigated 
and on the speed with which the mass matrices rotate, and these in turn 
depend on the theory or model one is considering. We examine below and in 
a companion paper |L4| some specific examples. 



4 Two Rotation Mechanisms 

In this section, we investigate the rotating fermion mass matrix in 2 spe- 
cific schemes at, as it were, opposite ends of the spectrum so as hopefully to 
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span most reasonable possibilities. At one end of the spectrum, we assume 
that there are no forces in nature other than those currently studied in the 
Standard Model so that all the rotation in the fermion mass matrix is driven 
by the electroweak Higgs terms via the empirical nondiagonal mixing (CKM 
or MNS) matrix. The generation phenomenon is thereby left unexplained. 
Fermion state vectors are to be specified by the FSD prescription (to con- 
trast with the scheme below which only works with the SSD prescription) 
in which the mass matrix is diagonalized at some chosen (low) scale. In 
the case of leptons, since some mass and mixing parameters are but poorly 
known, these specifications have yet to be supplemented by further assump- 
tions based on current popular views in a manner to be detailed later. We 
shall refer to this henceforth as the NSM scheme (for naively implemented 
Standard Model). At the other end of the spectrum, we consider a scheme 
we ourselves suggested |10| , which in turn was abstracted from our so-called 



Dualized Standard Model (DSM) || [12] although it need not strictly adhere 



to the duality concepts of that earlier work. In this scheme, generation it- 
self originates from a broken SU(3) gauge symmetry, the rotation is driven 
by Higgs bosons associated with the breaking of this symmetry, and even 
the nondiagonal mixing matrices as experimentally observed are themselves 
consequences of the mass matrix rotation. Here the SSD prescription for 
defining state vectors is essential. Given that in the DSM it is the rotation 
which gives rise to fermion mixing, and not the other way round as in the 
NSM case with the mixing driving the rotation, transmutation effects are 
usually, though not in all cases, larger in the DSM than in the NSM, and 
hence would be easier either to disprove or to confirm. 

Our first task is to evaluate the rotating mass matrix over a range of 
scales, say, up to the energy achievable in the foreseeable future of around 100 
TeV. For the NSM case, this in principle requires the solution of the equations 



(1.1) and ( |1.2|) together with any other renormalization group equations they 
are coupled to, which though standard is a little involved. For our present 
purpose of a first exploration, however, we may approximate by linearizing 
the equations and keeping only those terms actually driving the rotation, 
namely the term DD^U in (|1 . 1| ) and the term UWD in ( |1.2|) . We need 
as input from experiment the masses of the fermion states and the mixing 
matrices measured at low energies. Assuming that the scale in the FSD 
prescription chosen for defining the fermion state vectors coincides with the 
scale at which the empirical mixing matrices are measured, the rotating mass 
matrices can then be evaluated by iterating the two equations. For quarks, 
both the masses and the CKM matrix have now been determined to sufficient 
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accuracy for our purpose and pose no practical difficulty. This calculation has 
been performed but since it will not be of use in what follows in this paper, 
the result will not be presented. For leptons, however, the MNS matrix is 
still poorly known, while the Dirac masses of neutrinos on which their Higgs 
couplings depend are still almost entirely unconstrained. For these, therefore, 
we shall just insert for our exploration the most popular theoretical biases 



at present, namely for the MNS matrix the so-called bi-maximal mixing fl5 
version: 




U MNS = 1/2 1/2 l/y/2 , (4.1) 



and for the Dirac mass 7713 of the heaviest neutrino a value around the mass of 
the t quark. One shall need the Dirac mass m 2 of the second heaviest neutrino 
also, which is taken, for lack of any better choice, to be of the order of the 
charm mass. Since the equations have been linearized in approximation, 
the result of the rotation for charged leptons can be summarized simply by 
giving the rates of change of the off-diagonal elements of the ('hermitized') 
mass matrix with the logarithm of the energy, which are explicitly given as: 



dL 



dt 128tt 2 V 2 



( 2m\m e -^finl(m e + ra M ) ^m\{m e + m T ) \ 

\ ^75 m i( m e + m r) 7713(771^ + m T ) 2m\m r j 



(4.2) 

with V being the vev of the electroweak Higgs, i.e. 246 GeV. Putting in 
then the suggested values of the masses, one obtains that (fi\m\r) changes 
by about 5.5 x 10~ 3 , (e\m\r) by about 1.77 x 1CT 7 and (e\m\fi) by about 
1.06 x 1CT 8 GeV per decade change in energy. These results will be useful 
later for estimating the rates of various transmutation effects implied by the 
NSM scheme. One notices that in ( |4.2j ), (e|m|r) is proportional to m\ only by 
virtue of the assumed zero in the bi-maximal mixing matrix. If that element 
U e 3 of the MNS matrix is nonzero, then {e\m\r) would be proportional to 
7713, which in turn can lead on iteration to terms of order m\ in the element 
(e|m|/x), thus making both much larger than the above estimates, which 
represent therefore a sort of lower limit. 

For the DSM case, the calculation, though more intricate, has been more 



explicitly specified. As the scheme depends on only 3 parameters |L0| which 
have already been determined by fitting to some and then used to give very 
sensible predictions to others of the mass and mixing parameters of both 
quarks and leptons, no further information needs to be supplied. Indeed, 
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in making these fits and predictions, the rotating mass matrices for quarks 
and leptons had perforce already been evaluated, though not previously pre- 
sented. The result for the charged leptons which is of concern to us in this 
paper is now given in Figures ^. 

To check our result, without going through the details of the fit in JTD|, 
one needs only to iterate the following equation for a vector (x', y', z'): 



(4.3) 



with 

, x'(x' 2 - y' 2 ) x'(x' 2 -z' 2 ) 

X /Z _|_ ylZ _|_ Z IZ 

starting from some initial value (xj, yj, zj) with a coupling p, all given in eq. 
16 of |H|. For an accuracy of about a percent, some 500 iterations are needed 





per decade change in energy, where the vector (x',y',z') is to be normalized 
after every iteration. The resulting mass matrix is of the factorized form: 



x' 



rri = m T \ y' {x',y',z'), (4.5) 



z' 



where ttlt is a normalization depending on the fermion-type T, which may 
be taken as the mass of the heaviest generation state of that type. In Q4.5| ), 
the mass matrix is given with respect to a fixed (gauge) basis in generation 
space, but the matrix with respect to the physical basis as presented in, for 
example, Figure || is also easily evaluated with the (physical) state vectors 
determined in accordance with the SSD prescription outlined in section 2. 
Thus, the physical state vectors of the charged leptons in the gauge basis 
are, explicitly: 

v T = (0.9967,0.0760,0.0268) 

= (-0.0759,0.7741,0.6285) 
v e = (0.0271,-0.6285,0.7774). (4.6) 

For more details, the reader is referred to [|l(| and earlier references therein. 

One particular feature of the DSM scheme is that all fermions lie on 
the same trajectory of the rotating vector (x', y', z') so that at any scale the 
mass matrix for up and down fermion types have the same eigenvectors. This 
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means that the U and D matrices in ( |1.1| ) and (O) can be simultaneously 
diagonalized, causing thus no de-diagonalization to each other via the elec- 
troweak Higgs loop effect. We recall that it was the latter which was what 
drove the rotation in the NSM case. Here instead all the rotation comes from 
the Higgs bosons associated with the broken 'horizontal' symmetry. 

Comparing the actual values for the rotating mass matrices of the two 
schemes NSM and DSM, one finds the off-diagonal elements to be generally 
larger for the latter than the former, as we expected. Further, as can be 
seen in Figure |3|, the mass matrix in the DSM scheme has quite an intri- 
cate structure with, for example, the off-diagonal elements passing through 
zero whenever the energy scale equals the mass of a physical fermion state, 
meaning that the mass matrix is there diagonal. This is a consequence of 
the SSD prescription by construction, and would apply also to the NSM had 
one chosen there also the SSD instead of the FSD prescription. 



5 Examples of Transmutational Decays 

Once given the rotating mass matrix as a function of the energy scale, it would 
appear that one has in principle already enough information to evaluate 
any implied transmutational effect. However, because of the mass matrix 
rotation leading to, among other things, some rather unfamiliar kinematics, 
adaptations to usual calculational procedures are required, which are not, 
at first encounter, entirely straight-forward. We find that this is true even 
in evaluating, for example, the simple Feynman diagrams of Figure ||] for 
the photo-transmutation of leptons that we have calculated. Therefore, in 
order to avoid clouding immediately the essentially simple general picture 
of transmutation with technical details, we relegate detailed calculations to 



other papers, e.g. on photo-transmutation to a companion paper [fL4|| , and 
limit ourselves here to qualitative estimates of widths for some illustrative 
examples in transmutational decays, which being single particle effects are 
simpler to analyse theoretically. Also, to avoid the complications of quark 
confinement, we shall deal here with only lepton transmutations, which are 
already of considerable physical interest, being constrained by some very 
stringent bounds in experiment. 

(I) As a first example, consider the decays: 

Z° — U/*, (5-1) 
with l a ^ lp being different lepton states. As indicated already in Section 3, 
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the amplitudes of such off-diagonal transmutational processes are expected to 
be suppressed with respect to those of their diagonal counterparts Z° — > l a l a 
by a factor of the form (a\m\(3) /E, with E being a measure of the typical 
energy carried by the fermion line, which we can take here to be of the 
order of (half) the Z° mass Mz- The numerator (a\m\/3) represents the off- 
diagonal mass matrix element linking the lepton states l a and lp at the scale 
of the reaction, which is again Mz in the present case. Given any choice 
for l a and lp, (a\m\j3) can be read off for the DSM in Figure ||. Thus, for 
example, for the decay Z° — ► t~ jX^ the relevant matrix element (/x|m|r) at 
the scale Mz = 91GeV is seen to have value about —0.11 GeV. For the NSM, 
the element (fj,\m\r) was said in Section 3 to vary by about 0.5 percent per 
decade change in energy. According to the FSD prescription adopted there, 
we have to first fix a scale at which the mass matrix is diagonal. If this is 
taken, for lack of any obvious better choice, to be at the r mass, we obtain 
for (/i|m|r) at Mz a value of about 0.01 GeV. To estimate the width for the 
decay Z° — > r~/i + , the easiest way would be to compare it with either of 
the diagonal decays Z° — > t~t + or Z° — > fi~fi + . Given that the masses of 
the final state leptons are so much smaller than the decaying Z° mass, one 
can neglect the small differences in kinematics between the different decay 
modes and write simply: 

r(z° — > rVj K/iHr)! 2 

r(Z o _^ r - T +) ~ M | ■ 

Putting in the above estimates of the matrix elements and the empirical 
value 3.36 percent jTJI for the branching ratio of the diagonal r~r + mode, 
one obtains then for the branching ratio of the transmutational mode t~/j, + 
an estimate of about 4 x 10" 8 for DSM and 4 x 10" 10 for NSM. Both of these 
estimates are much below the present empirical bound of 1.2 x 10~ 5 listed 
I6|| . Similarly, estimates for the branching ratios of the other modes are 



m 



obtained, giving for the r~e + mode a branching ratio of 4 x 19~ 9 and for 
the n~e + mode 1.6 x 10 -11 for DSM. Again these are way below the present 
empirical limits listed in JTH] of respectively 9.8 x 10~ 6 and 1.7 x 10~ 6 . The 



estimates in the NSM for these last two modes are minuscule due to the very 
small values for the mass matrix elements linking e to fi and e to r, as will 
also be true for the following examples. They shall therefore henceforth be 
ignored except, for a particular reason to be explained, for (IV) below. As 
one can easily see, the reason that the transmutation rates are all so small 
in Z° decay is that the off-diagonal effects are much suppressed by the high 
energies involved which is here of the order of the Z° mass. 
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(II) As a second example then, let us go to the other extreme in energy 
and consider the decay: 

7T° — ► /_Te + (5.3) 

which in 7r°-decay is the only kinematically accessible transmutational lep- 
tonic mode. The same arguments as before suggests that we compare this 
with the diagonal mode 7r° — > e~e + giving a suppression in the amplitude 
of the former with respect to the latter by a factor of the form (e\m\fi) /m^ 
where the off-diagonal mass matrix element is to be evaulated at the scale 
77V Here, the spins, parities and masses of the particles involved being such, 
a correction for the difference in kinematics between the two decays is war- 
ranted, which is made most easily by comparison with the analogous charged 
Ti decays, giving thus: 

rfrQ — > ^- e +) \(e\m\^T(n±^^v,) 

r(7r° — > e-e+) ~ 77i2 r(vr± — > e±v e ) ' 1 ' ' 

The matrix element (e|m|/x) at scale m n o = 135 MeV obtained from the 
calculation giving Figure || is about 3 x 1CT 4 GeV. Hence, putting in the 
empirical value of 7.5 x 10~ 8 for the branching ratio of the diagonal mode 
7T° — > e~e + and that of 1.23 x 10~ 4 for the branching ratio of r(7r ± — > e ± i/ e ) 



16| , we have an estimate for the branching ratio of the transmuational mode 
7t° — > ^Te + of about 2.9 x 1CT 9 . This is barely an order of magnitude 
below the experimental bound of 1.72 x 10~ 8 [[UJ. Indeed, one is saved here 
from a violation of the experimental limit only by the exceptionally small 
value of the off-diagonal mass matrix element due to the proximity of the 
reaction scale m^o to the muon mass where the mass matrix is constrained 
to be diagonal by the SSD prescription. 

(Ill) As a third example, let us consider at an intermediate scale: 

ip — ► Ijp. (5.5) 

Following the same procedure as above, one obtains: 

rty— >t-a* + ) |(/i|m|r)| 2 



(5.6) 



The matrix element (/x|m|r) for DSM as read from Figure ^ at the scale m^> 
is about 0.03 GeV. Combined with the empirical value of about 6 percent 



[ 16fl for the branching ratio of ip — y H one obtains an estimate of about 
6 x 10~ 6 for the transmutational r^/j + mode, which is larger than in either of 
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the two earlier examples. The corresponding estimate for the NSM is about 
1 x 1CT 8 . Although no empirical limit is given in for this particular mode, 
the estimates seem close enough to the sensitivity of present experiment to 
make a search for it worthwhile. Similar arguments applied to the modes 
if) — > r~e + and ip — ► fi~e + yield the estimated branching ratios 1.6 x 10 -7 
and 6 x 10 -11 . For these modes again, no empirical limits are given in (JTSj but 



the estimates seem in any case beyond present experimental sensitivity. The 
same anlysis can also be applied to T decay giving for the branching ratios 
of the modes r~fi + ,r~e + , and /i~e + respectively the following estimates: 
1.6 x 10~ 6 , 6 x 10~ 8 , and 7 x 10 -11 . Again, no empirical bounds for these 



modes are given in IfTfi] , but some, especially the r/x mode, will be worth 
searching for in future in e.g. 5-factories. 

(IV) Lastly, consider transmutational fermion decay, of which the prime 
examples would be: 

fT — > e~7 (5.7) 

and 

jjT — > e~e + e~ (5-8) 
proceeding via the Feynman diagrams in Figure 01. The bound as given in 





(b) 



Figure 4: Transmutational \i 



e 7 and \x~ 



e e + e decays 



| 16| for the branching ratio of the decay (|5.7| ) is an impressive 4.9 x 1CT 11 , 
and of the decay (|5.8| ) is 1 x 10 -12 . And these coming on top of a main 
decay mode of /i, namely \x — >■ evi>, proceeding by weak interactions with a 
width of only around 3 x 10 -10 eV, are a very stringent constraint indeed on 
any scheme which violates muon number. It is thus gratifying for us to note 
that transmutation in the DSM scheme survives this very dangerous-looking 
hurdle automatically. This comes about because the Feynman diagrams H 
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for the transmutational decays ( p.7| ) and ( |5.8| ) are to be evaluated at the 
scale of the \i mass, and at that scale the mass matrix in the DSM scheme is 
diagonal. That it is so by virtue of the SSD prescription for defining lepton 
states was noted already in preceding sections and can be seen explicitly in 
Figure |3| by the vanishing of the off-diagonal mass matrix element (e\m\n) 
at the \i mass. Hence the DSM scheme actually predicts zero probability for 
this transmutational decay. 

The same observation applies also to transmutational decays of r, for 
example: 

t~ — > e"7, /x"7, (5.9) 

r~ — > e~e + e~ , e~/i + /i~, /i~e + e~, (5.10) 

since for these decays the analogous Feynman diagrams are to be evaluated 
at the scale of the r mass, and at that scale the mass matrix in the DSM 
scheme is again diagonal giving thus zero transitions. This is also fortunate, 
for although not quite as stringent as those on ( |5.7|) and (|5.S| ), the bounds 
on the branching ratios of (|5.9|) and (|5.10 ) as given in [16] are still only of 



order 10 -6 for a total r width of around 3 x 10~ 3 eV, which would be far from 
trivial to satisfy otherwise. 

To illustrate how sensitive a test these lepton decays are to transmutation 
models, let us examine them in the NSM scenario. Here, as noted before, the 
rotation of the mass matrix is slower than in the DSM case giving thus usually 
smaller transmutation rates, especially for transmutations between \x and e. 
However, the mass matrix is required by the FSD prescription adopted to be 
diagonal only at some fixed mass scale. Suppose we choose this scale as the r 
mass so as to make the transmutation rates vanish for both the decays ( |5.9| ) 
and ( 5.10Q and so satisfy the experimental bounds on them automatically. As 



the energy lowers to the scale of the \i mass, however, the mass matrix will 
have rotated and acquired nonzero off-diagonal elements. Indeed, according 
to (|4.2| ) the element (e|m|yu) would have acquired on running from the r 
mass a value of around 1.3 x 10 -8 GeV, giving thus possibly nonzero rates 
to the decays (|5.7| ) and (|5.8|) . To estimates these rates following the previous 



procedure, let us compare for example ( |5.8| ) to the known main decay mode 
of /i, namely: 

pT — > e"z/ e z/ M , (5.11) 

as given to first order by the Feynman diagram of Figure |5]. The mass of 
the electron being so much smaller than that of the muon, we can ignore the 
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Figure 5: Feynman diagram for the decay fi — ► e v^v^ 



small difference in phase space between the two decays and write: 



r(/i~ — >e~e + e~) \(e\m\fi)\ 2 e 2 /ml 
r(/i" — > e-u e u^) ml g$/M& 



(5.12) 



where we have taken for a typical momentum for the decay to be of order of 
the muon mass. Inserting the above NSM value for (e\m\/j,) and the empirical 
values for the other quantities appearing in (|5.12 ), one obtains a value of 



about 1x10 3 for the ratio. This would be disastrous, being some 9 orders too 
large compared with the experimental bound given in [16| of order 1 x 10 -12 . 



That this ratio turns out to be so large despite the very small NSM value for 
(e|m|/x) is because of the large W mass occuring in the diagram for the normal 
decay mode of /x, or in other words, because we are in ( |5.12| ) comparing what 



is basically, in spite of being off-diagonal, still a photon-mediated process to 
a weak decay. The above blatant violation of the experimental bound is of 
course not the fault of the Standard Model but only of the naive fashion that 
it has been implemented in the NSM scheme since it could be easily avoided 
by, among other means, choosing a smaller value for the Dirac mass m 2 of 
the second heaviest neutrino which is still at present entirely unconstrained 
by experiment. Nevertheless, this example shows how sensitive these lepton 
decays £1X6 clS £L test of transmutation models, and how delicately the DSM 
scheme has managed to survive it. 

Obviously, what really mattered for passing that test was what we called 
the SSD prescription for defining the lepton state vectors which stipulates 
that the mass matrix is to be diagonal at the mass of every lepton state. 
Had one chosen in the NSM scheme to define fermion states in this way 
instead of the FSD prescription adopted, it would have passed the test also. 
With the FSD prescription, however, it would not have mattered if one had 
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chosen the scale for diagonalizing the mass matrix to be m M instead of m T , 
for although one would then pass the test for the decay //~ — > e~e + e~, 
one would fail it for such decays as r _ — > \T . Indeed, running now 
from the [i mass, one would obtain a non-diagonal mass matrix at the r 
mass with off-diagonal elements of the same order. The empirical bounds on 
transmutational r decays are less stringent, e.g. the present bound on the 



branching ratio of t~ — > is only 1.9 x 10 ||16|| , but that would 

still mean a violation of the bound by some 6 orders in magnitude by the 
above estimate. One could of course impose in the NSM a further condition 
on the FSD prescription stating that the mass matrix should stop rotating 
altogether once below the chosen diagonalization scale, but one does not 
know of a good theoretical reason for doing so. 

Firm conclusions cannot be drawn until detailed calculations for the var- 
ious decay widths have been performed, which to us seem quite feasible in 
the near future. Tentatively, however, one can claim that the DSM pre- 
dictions for transmutation survive present experimental bounds in all cases 
so far studied, while the NSM (with FSD) survives except perhaps for the 
noted /i or r transmutational decays. In most cases, the predicted rates are 
some way below present experimental limits, but for the decays ir° — > e ± fi T , 
if) — > / u ± r =F and T — > / u ± r =F the estimated rates may soon be, if they are 
not already, within range of experimental sensitivity. 



6 Remarks 

The importance of possible fermion flavour violation has long been recognized 
and subjected to rigorous test by such experiments with ultra-high sensitivity 
as fi — > e7 decay. The negative outcome to-date of all these tests conveys 
the impression that such violations can exist, if at all, only to a very small 
extent. However, that conclusion was reached in the days before it was 
realized that the fermion mass matrix actually rotates with the energy scale 
so that the fact that flavour is conserved at one scale, or that the mass matrix 
is there diagonal in the flavour states, does not necessarily mean that it will 
remain so at some other scales. That being the case, the accuracy, though 
impressive, of such experiments as /i — ► searches, is not enough to rule 
out flavour violation in general but has to be supplemented by experiments at 
other energies. As our analysis above has shown, under some circumstances 
the mass matrix can acquire quite appreciable off-diagonal elements without 
being noticed in existing experiments. Hence, regardless of which scheme for 
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mass matrix rotation one favours, it would be worthwhile for experimenters 
to perform routine checks under varying circumstances for the diagonality of 
fermion mass matrices by looking for effects of the type labelled in this paper 
as transmutations. With a rotating mass matrix, diagonality is no longer as 
sacrosanct as once believed. 

As argued in the Introduction, that the mass matrix rotates forces on us 
a refinement of our definition of fermion flavour states. Given the sensitivity 
of such experiments as /i~ — > e~e + e~ searches, even the small rotations 
which may occur between the mass scales of two successive generations is 
enough to give huge violations of the experimental limits, as we have shown 
in the preceding section. It seems thus necessary either to define the fermion 
states in such a way as to ensure that the mass matrix is diagonal at the 
mass scales of all the fermion states, as was done in the SSD prescription, or 
else in a prescription like FSD to suppress altogether the rotations between 
fermion states, if not by some theoretical mechanism then by decree. For us, 
this conclusion has come as an agreeable surprise in that the SSD prescrip- 
tion which was originally invented just for internal consistency of our DSM 
scheme should turn out now to have the unexpected virtue of guaranteeing 
consistency with experiment on \i — > and (i~ — > e~e + e~ decays, which 
would otherwise be a difficult hurdle to survive. As we shall see in iTl 



the same property of the SSD prescription will save one also from certain 
unwanted pole structures in the amplitudes for the photo-transmutation of 
leptons. 

The rotation of its mass matrix being weaker, the NSM scheme will give 
generally smaller transmutation effects than the DSM, i.e. provided that it 
has been amended with something like the SSD prescription for defining 
fermion states to avoid the noted probable inconsistencies. Indeed, in all 
the examples so far studied NSM transmutation effects seem to be beyond 
present experimental sensitivity. That being the case, if in the routine check 
on transmutation bounds an experiment comes up with an effect, then it is 
probably due to a rotation of the mass matrix driven by some forces of nature 
other than those currently considered in the Standard Model. This would be 
exciting in any case, besides possibly opening up a new window for probing 
into the generation mystery, as we have already said in the Introduction. 

For us advocates of the Dualized Standard Model (DSM), the results of 
this paper and its companion []H are something of a relief, at least temporar- 
ily. Since the whole DSM scheme relies heavily on the concept of a rotating 
mass matrix driven by dual colour forces, on which all its seemingly sensi- 
ble predictions on the mass and mixing patterns of fermions also depend, it 
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would be a little disastrous if the relatively fast mass matrix rotation required 
for the scheme's success predicts at the same time much larger transmutation 
effects than can be accommodated by present experiment. And, in contrast 
to flavour- changing neutral current (FCNC) effects |T7], [H| [TB[ which also 
pose some hurdles for the scheme to overcome but depend on an unknown 
mass parameter, there is in transmutations no free parameter which could 
be adjusted as a loophole for escape. It is thus fortunate that a disaster 
does not seem so far to have happened although in some cases it has come 
rather dangerously close. Instead, one obtains as a bonus a whole new class 
of phenonmena, with quantitative predictions which, while remaining below 
present bounds, are still not that far off as to be entirely inaccessible to exper- 
iment. However, much more work will still be needed first, to make sure that 
disasters will not happen in circumstances not yet examined, and secondly, 
to identify specific cases where the scheme's predictions can be checked by 
experiment, if not immediately, at least in the not too distant future. 

In general terms, what we have done in this paper is to raise some basic 
questions of concept forced on us by the rotation of the fermion mass matrix 
and attempt to answer some of them. Although in answering these questions 
we may have raised as many more, we hope at least to have clarified a little 
some of their implications. 
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